Superfiuid to insulator phase transition in a unitary Fermi gas 
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We study the evolution of the energy gap in a unitary Fermi gas as a function of temperature. 
To this end we approximate the Fermi gas by the Hubbard lattice Hamiltonian and solve using the 
dynamical mean-field approximation. We have found that below the critical temperature, T c , the 
system is a superfiuid and the energy gap is decreasing monotonously. For temperatures above T c 
the system is an insulator and the corresponding energy gap is monotonously increasing. 
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Introduction - Dilute Fermi gas, characterized with in- 
terparticle distance much larger than the effective range 
but much smaller than the scattering length, |o a | 3> 

\J r eff, has been the subject of intense theoretical 

and experimental research in the last few years [l[ . The 
interest in this system stems from its universal properties 
that become independent of its actual constituents as the 
scattering length diverges at unitarity (\a s \ — ► oo) and 
depend only on the particle density. In the weak cou- 
pling regime (a s small and negative) the ground state of 
a Fermi gas is a BCS superfiuid. In the strong coupling 
limit (a s small and positive) the fermions are bound in 
pairs that form a Bose-Einstein condensate (BEC) for 
temperatures below the critical temperature. At unitar- 
ity the system is in between these two limits and exhibits 
a distinct behavior which can be classified as a new type 
of superfluidity, characterized by an admixture of bosonic 
and fermionic features ||. 

The phenomenon of superfluidity in Fermi systems is 
associated with the occurrence of off-diagonal long range 
order, A = U{Tc^(0 + )ci(0)), and the existence of a 
gap Ag a p in the single particle excitation spectrum. In 
general, the order parameter An and the gap A gap are 
independent quantities. However, for weakly interacting 
fermions, in the BCS regime, one finds A gap = Ao. 

In an intriguing paper, Bulgac et. al. Q used quantum 
Monte Carlo technique to study the evolution of A gQp as 
a function of temperature assuming a quasi-particle spec- 
trum. They have found that in contrast with the BCS 
theory, where A gap vanishes at the critical temperature 
T c , for unitary Fermi gas the magnitude of A sap (T c ) is 
about two-thirds of the zero temperature gap. Even more 
striking is the fact that at about T c the gap's derivative 
flips sign, i.e. A gQp grows for T > T c . 

These results call for a better understanding of the ex- 
citation spectrum of the finite temperature unitary Fermi 
gas. The aim of this work is to study these aspects 
of the system using the dynamic mean field approxima- 
tion (DMFA) (3, If. In the DMFA, a lattice problem is 



mapped into a self-consistent embedded impurity prob- 
lem. In the limit of infinite spatial dimensions d — > oo 
this mapping becomes exact due to the localization of the 
self-energy Q. For 3D fermions which we consider here, 
DMFA can be regarded as a simplification in which a 
purely local self-energy is assumed, S(fe,io; n ) m T,(iuj n ) 
(hat denotes a spinor matrix). The validity of this as- 
sumption for unitary Fermi gas has been examined in 
0, [I] , where the problem was approximated by the lat- 
tice Hubbard Hamiltonian, 

a riri' n 

and the continuum limit was realized by reducing the 
lattice filling to zero. It was found @ that the DMFA 
results agree remarkably well with those of full quantum 
Monte-Carlo simulations (QMC) [1 M, M, El , yielding 
£ w 0.44 for the ratio between the energy per particle 
of the interacting and free systems (E/N = ^Epc), and 
A w QMEp. In this letter we use the DMFA to study 
the finite temperature excitation spectrum of the unitary 
Fermi gas. 

DMFA - Using the Nambu formalism, the DMFA 
single-site impurity effective action takes the form 

Seff = ~ f dT ( dT'&{T)Qv l {T-T')*{T') 
JO Jo 







U / dTc\{T)c^T)c\{T) Cl {T) 



(2) 



where (3 = l/T is the inverse temperature, = 
(cj,C|) are the Nambu spinors, and the bath's Green's 

function Qq is determined through the self-consistency 
condition that the impurity Green's function Q{t) = 
— (T3 r i (r)^(0))s ef/ coincides with the site-diagonal lat- 
tice Green's function calculated with the self-energy 

E(»W n ) = G^iWn) - G^ii^n) ■ 

We use the direct diagonalization method of Caffarel 
and Krauth [l3| to solve the DMFA. In this approach the 
impurity action is mapped into the Anderson Hamilto- 
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where the interaction of the fermionic field c a with the 
auxiliary bath fermions a; CT generate Qq. This goal is 
achieved by choosing the parameters of the Anderson 
model ei,Vi,Di to minimize the difference between the Go 
and G^ nd . In this work we use 4 — 5 auxiliary fermionic 
fields. For lattice filling n = 0.1 which we consider here, 
taking this number of auxiliary fields yields an accuracy 
of about 1% for the thermodynamic quantities 0]. 

The Excitation Spectrum - The determination of real- 
frequency quantities such as the spectral function or the 
excitation spectrum faces severe limitations in QMC sim- 
ulations where only imaginary time/frequency data are 
obtained directly. Trying to overcome this limitation 
Bulgac et. al. [3| have calculated, using a QMC sim- 
ulation, the susceptibility 



x(fc) 



f drG{k,T) = -\y —G{k,iLu n ) , (4) 
Jo P ' tu n 



where G{k 1 r) is the Green's function and to n — 

(2n + l)Tr/f3 are the Matsubara frequencies. For an 

independent- (quasi) particle spectrum the response ((4]) 
can be easily evaluated, 
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where E k are the single-(quasi)particle excitation ener- 
gies. Exploiting this observation, they have fitted the cal- 
culated susceptibility to the formula ([5|) assuming, given 
the chemical potential /x, the spectrum 



\J (a qp e k + T, qp - fj,) 2 



A 2 

qp 



(6) 



treating a qp ,Y> qp , A qp as free parameters. These param- 
eters stand for the effective mass m* = m/a, mean field 
potential ^ qp , and the "pairing" gap A gap — A qp . e k is 
the free-particle kinetic energy. 

Solving the DMFA equations with the direct diagonal- 
ization method [l3| the spectral function can be obtained 
directly from the impurity model but in the form of a set 
of delta functions. Since we are limited to a finite and 
rather small number of orbitals in the effective bath it is 
difficult to extract quantitative information from it. Con- 
sequently one has to adopt a different strategy in order 
to study the real frequency properties of the system. 

Consider the occupation probability 



/(fc) = G(M^ 



G{k, icu ri 



(7) 



the Green's function derivative at r = H 



C(fe) = 



dG(k,r) 



T=0+ 
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e^" ' ~ iu n G(k,iuo n ) , 



and the susceptibility x(&) defined above ((4]). For an 
independent-(quasi)particle Green's function, 



G qp (k, iuj n ) 



(9) 



(iuj n - E k )(iuj n +E k ) ' 
these quantities can be manipulated to yield the relation 



E k = \ — 
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This procedure is a generalization of 0] that avoids, 
however, the need to invert Eq. ([5]). In the DMFA 
x{k), f(k)X(k) can be easily calculated through the 
Matsubara sums (|3]), ([7]), and ([H]). Once we have per- 
formed these sums the value of E k can be evaluated (| 10[) , 
regardless of the original assumption about the nature of 
the excitation spectrum. Strictly speaking, only for a 
limited number of cases the identification of Eq. (JTOJ) 
with the quasi-particle excitation spectra is exact. Nev- 
ertheless, in the following we shall refer to it as the quasi- 
particle excitation spectrum. 




FIG. 1: (Color online) The quasi-particle spectrum as calcu- 
lated from Eq. |TlO)) for a unitary Fermi gas at T — 0.38-Ef 
(T c ~ 0.16-Ef) and lattice filling n = 0.1. The minimum in 
the graph corresponds to the quasi-particle gap A qp . 

In Fig. Q] the quasi-particle energy (fTTTj) is plotted as a 
function of the free-particle kinetic energy e k for lattice 
filling n = 0.1, at temperature T — 0.38Ep, beyond the 
phase transition temperature which in our calculation is 
T c « O.WEp. From the figure it can be seen that E k 
exhibits an excitation spectrum typical for a gapped sys- 
tem, to which we shall refer as an insulator. The gap can 
be evaluated directly from the graph as the minimum 
of E k . Fitting the empirical formula © one can repro- 
duce the excitation spectrum (|10[) very accurately. Using 
this procedure we get an estimate for A qp even when the 
minimum of (|10[) is outside the band and we also get an 
estimate for the effective mass. In Fig. [5] we plot for a 
unitary Fermi gas the gap A qp and m/m* as a function 
of T at lattice filling n = 0.1. Also plotted are the energy 
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per particle E, the chemical potential fi, and the order 
parameter A . The phase transition from a superfluid 
to a normal phase is associated with the vanishing of the 
order parameter and a jump in the heat capacity. Such 
transition is easily located in Fig. [2] at T » O.WEp. Fol- 
lowing the evolution of A qp with T, we see that at low 
temperatures A qp is a decreasing function of tempera- 
tures up to T c . At which an abrupt change is observed 
and for T > T c the gap is increasing with T. This dis- 
continuity in the derivative dA qp /dT is a clear indication 
that the quasi-particle gap has a different meaning in the 
two phases. For T < T c it can be associated with the su- 
perfluid gap, Ag a p. This interpretation, however, is lost 
for T > T c . 
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FIG. 2: (Color online) The phase transition from a superfluid 
to a normal phase in a unitary Fermi gas. The energy per 
particle E is shown by squares, the order parameter Ao by 
circles, the chemical potential \i by up triangles, the quasi- 
particle gap Aqp by down triangles, and the bare to effective 
mass ratio m/m* by diamonds. 

A better understanding of the quasi-particle gap at 
T > T c can be achieved by frustrating the superfluid 
phase at T < T c . Within the DMFA this goal can be 
easily achieved by forcing particle number conservation 
in the effective impurity action, i.e. by setting Di = 
in ([3]). In Fig. [3] we present the results of such calcula- 
tion for the quasi-particle gap. At very low temperatures 
T < T, 



pairing 



0.02Ep the frustrated solution exhibits 
a gap of about 0.13Ep, which corresponds to the coex- 
istence region of a metal phase and a pairing (insulator) 
phase. This phase transition in the normal, unstable, 
phase around the unitarity limit was identified at T = 
by Keller et. al. 1 141) for an infinite dimensional system. 
Toschi et. al. [15| have studied the finite temperature 
phase diagram of the frustrated solution, and found that 
although there is a smooth transition between the metal 
and insulator phases at temperature above the metal- 
pairing critical point (which is much lower than T c ) the 
properties of the system depend strongly on the strength 
of the coupling constant. Above T pa i r i ng we see a drop in 
the normal phase gap which then grows almost linearly 



with T until T c . It is interesting to note that the de- 
creasing superfluid gap and the increasing normal phase 
gap coincide at T c . It is evident that for T > T c the 
thermodynamic stable gap A qp shifts from describing the 
superfluid solution into the normal phase solution. From 
the figure it seems that the gap saturates at higher tem- 
peratures. Being limited by the band width there is no 
point in carrying our calculations to higher values of T . 
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FIG. 3: (Color online) The quasi-particle gap as a function of 
T. The thermodynamic stable solution is shown by triangles. 
The frustrated superfluid solution (the normal solution) is 
presented by diamonds. 

We can achieve further insight into the behavior of the 
unitary Fermi gas by inspecting the self-energy, E(iw n ), 
at T > T c , see Fig. It can be seen that the real 
part of S(zw„) is essentially constant. The imaginary 
part decrease asymptotically as l/u> n but tends towards 
a finite value as u> n — ► 0. Consequently S(u/ n ) can be 
roughly represented in the form, 



S 



liO n + i&O 



(11) 



where 770 characterize the asymptotic behavior of ImS 
and #0 is used to model the low frequency behavior. The 
dashed line in Fig. 0] is a fit of (fTT]) to the calculated 
self-energy at T = 0.365Ep- The best fit parameters are 
So = -0.20.Ef, 77 = 0.82E F , and 6 = 1.04E F . Using 
(fTTj) we get an analytic model for the thermal Green's 
function. The upper plane part of this Green's function 
can be related to the retarded Green's function yielding, 

G r (lo, k) = ' 



fi - efe - So - 



(12) 



U+10Q 



In this model the Green's function contains two poles 



- iOo ± V( e k +i0o) 2 +^Vq 



where ey~ = + Sq — /i, and can be written as 
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FIG. 4: (Color online) The self-energy of a unitary Fermi gas 
at T = 0.365-Bf and lattice filling n = 0.1. The real part 
of E and imaginary part multiplied by ui n are plotted with 
thick lines. The dashed line is the simplified function (|lll) 
with best fit parameters £o = — 0.20-Ep, rjo = 0.82Ep, and 
0o = 1.04.Ef. 



The two poles come close to each other as approach 
the Fermi surface. At e*. = [i — T,q we g et a simp le ex- 
pression for the gap, Auj = u> + — lu_ = \J Arj^ — 9^ . This 
expression implies that for 770 > $o/2 there is a real gap 
in the spectrum and the system can be characterized as 
an insulator whereas for 770 < 9q/2 there is no gap and 
we may characterize the system as a Fermi liquid. For 



agreement with the quasi-particle gap A qp rj O.QEp, see 
Fig. [3l Inspecting the residues R± it is clearly seen that 
as ej. pass trough the Fermi surface the power is shifted 
from w_ to u>+, 
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the example in Fig. |4j A 
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This simple, but rather exact, parameterization implies 
the existence of a gap in the excitation spectrum and 
supports the conclusions we drew from the quasi-particle 
spectrum (fTU)) . 

Conclusions - Using the dynamic mean field approxi- 
mation we have studied for a unitary Fermi gas the evo- 
lution of the quasi-particle gap with temperature. We 
have found, in accordance with QMC calculations 0], 
that in the superfluid phase the gap decreases up to T c 
and then starts to rise. We have shown that at T c there 
is a sharp transition in the gap's slope. This transition is 
associated with a shift from the superfluid to the normal 
phase gap. We have demonstrated, by frustrating the 
superfluid solution, that the normal phase insulator gap 
is much smaller than the superfluid gap at low tempera- 
tures. The insulator gap grows with increasing tempera- 
ture and the two gaps coincide just at T c . The connection 
between the lost of coherence and the gap crossing is not 
yet clear. 
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